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Abstract 

We investigate the neutral Higgs sector in the next-to-minimal supersymmetric standard 
model (NMSSM) with explicit CP violation at the one-loop level by using the effective 
potential method. In general, explicit CP violation is possible at the tree-level in the Higgs 
potential of the NMSSM, which may possess a complex phase. The tree-level Higgs potential 
can be made CP conserving by assuming that all the relevant parameters are real. However, 
the CP-conserving Higgs potential at the tree level may still develop complex phases at 
the one-loop level through radiative corrections. These complex phases exhibit explicit CP 
violation through the mixings between the scalar and pseudoscalar Higgs bosons. One of 
the complex phases arise in the neutralino sector. For a particular choice of the relevant 
parameter values, the scalar-pseudoscalar mixings are found to strongly dependent on this 
phase. Meanwhile, for a wide region in the parameter space, the ratios of CP mixing elements 
of the neutral Higgs boson mass matrix to CP non-mixing ones increase weakly against this 
phase. Also, the neutral Higgs boson masses are roughly stable against the variation of this 
phase. 



I. Introduction 



Supersymmetry (SUSY) must be broken in order to be realized in nature. In various supersym- 
metric models, the SUSY breaking may be accomplished by introducing soft SUSY breaking 
terms. The introduction of soft SUSY breaking terms may also invoke irremovable complex 
phases. These complex phases may contribute to the mixing between the scalar and pseu- 
doscalar Higgs bosons, thus account for explicit CP violation. In the minimal super symmetric 
standard model (MSSM), there can be as many as eleven non-trivial phases associated with 
the scalar fermions, the charginos, and the neutralinos. However, at the tree level, the MSSM 
can avoid explicit CP violation by adjusting the complex phases away via redefinition of various 
participating fields in the model. It has been found that, at the one-loop level the MSSM can ac- 
commodate explicit CP violation in its neutral Higgs sector [1]. Many authors have investigated 
the explicit CP violation in the MSSM in the radiatively corrected Higgs potential [1-4], where 
five complex phases (three from the scalar fermion sector of the third generation, one from the 
charginos, and one from the neutralinos) can participate in the scalar-pseudoscalar mixing. 

Unlike in the case of the MSSM, explicit CP violation can be realized in the next-to-minimal 
supersymmetric standard model (NMSSM) even at the tree level. At the tree level in the Higgs 
potential of the NMSSM, a nontrivial complex phase emerges after redefining the three Higgs 
fields [5]. This complex phase from the tree-level Higgs potential may persist in the one-loop 
effective potential if the scalar top quarks are degenerate in mass. In Ref.[5], assuming the 
degeneracy of the scalar top quark masses, the scalar pseudoscalar mixing is investigated in the 
NMSSM with explicit CP violation. It has been shown that large mixings between the scalar 
and pseudoscalar Higgs bosons may be realized as the vacuum expectation value (VEV) of the 
neutral Higgs singlet approaches to the electroweak scale [5]. 

We have elsewhere investigated the implications of the mass splitting between the scalar 
quarks of the third generation in the NMSSM at the one-loop level with explicit CP violation 
[6]. In this case, two additional complex phase are induced from the scalar top and scalar 
bottom quark masses. We have then extended our investigations to the chargino sector in the 
NMSSM with explicit CP violation [7]. An additional complex phase other than the tree-level 
phase appearing from the mixing between the chargino masses is shown to contribute the scalar- 
pseudoscalar mixing, too, at the one-loop level. 

In this paper, we investigate the possibility where the tree-level Higgs potential is CP conserv- 
ing while explicit CP violation occurs in the Higgs sector of the NMSSM through the one-loop 
radiative corrections. We assume that the tree-level complex phase is zero. All the parameters in 
the NMSSM Higgs potential at the tree level are assumed to be real. Then, the possible sources 
of complex phases contributing to the scalar-pseudoscalar mixing are the radiative corrections 
from the scalar top, the scalar bottom, the scalar tau lepton, the chargino, and the neutralino 
sectors. In fact, in the one-loop effective potential for each sector may be embedded a complex 
phase. We investigate the effects of these five complex phases in the radiatively corrected Higgs 
potential of the NMSSM. One of the complex phases comes from the neutralino sector. We pay 
our attention on this phase. The dependence on this phase of the scalar-pseudoscalar mixings 
are investigated, while the values of other parameters are fixed. The scalar-pseudoscalar mixings 
are strongly affected by this phase, at least for a specific set of parameter values. 

Our paper is organized as follows. In Sec. II we briefly describe the Higgs sector of the 
NMSSM at the tree level. In Sec. HI, for the scenario of explicit CP violation, radiative 
corrections from the scalar top quark sector, the scalar bottom quark sector, the scalar tau 
lepton sector, the chargino sector, and the neutralino sector are explicitly calculated using the 
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effective potential method. Numerical analysis and results for the explicit CP violation at 
the one-loop level are presented in Sec. IV. Conclusions and some remarks are given in Sec. 
V. Appendices supply the formulae for the elements of the mass matrix of the neutral Higgs 
bosons, which contain the radiative corrections from each sector of the scalar top quarks, the 
scalar bottom quarks, the scalar tau leptons, the charginos, and the neutralinos. 

II. The Higgs sector 

The NMSSM has two SU(2) doublet superfields Hi = {H^,H-) and H2 = {H+,H^) with 
hypercharges —1/2 and 1/2, respectively, and a neutral SU(2) singlet supcrficld N with zcro- 
hypercharge. For the fermion matter fields, we take into account only the third generation of 
quarks and leptons. The general form of the NMSSM superpotential may be expressed as 

W = htQH2t% + hbQHib% + KLHiT'k + XNH1H2 - ^iV^ , (1) 

where H1H2 = H1H2 — H^ H2 and all the coupling coefficients arc dinicnsionless. The chiral 
superfields Q and L contain the left-handed quark and Icpton doublets, and t^j, 6jj and r]^ are 
the charge conjugate of the right-handed t quark, h quark, and r lepton, respectively, and the 
parameters ht, hh, and hr are their respective Yukawa coupling coefficients. 

In the NMSSM [8, 9], the tree-level Higgs potential may be decomposed into three parts as 

V*^ = Vd + Vf + Vs , 

where 

Vn = f{HldHi+HlaH2f + ^f{\H2f-\Hif)\ 

Vf = |A|2[(|i?i|2 + \H2\^)\N\^ + \HiH2\'^] + \k\^\N\^ - iXk*HiH2N*^ + H.c.) , 

Vs = ml^\Hi\^ + ml^\H2\^ + m%\N\^-{XAxHiH2N + hAkN'^ + R.c.) , (2) 

with gi and g2 being the U(l) and SU(2) gauge coupling constants, respectively, and a = 
{a^,a^,a^) are the Pauh matrices. The soft SUSY breaking Vs has two additional parameters 
Ax and A/^, both of mass dimension, and three soft masses m//-^, m^^^, and mjv- Generally, A, 
k, Ax, and Aj. can be complex. Among them, A^a and kA^ can be adjusted to be real and 
positive by redefining the phases of H1H2 and N . Thus the tree-level Higgs potential may have 
at most only one physical phase. It may be chosen to be in Xk* = Afce**^. This phase may 
allow for explicit CP violation in the NMSSM. It has been found that the scalar-pseudoscalar 
mixing in tree-level Higgs potential of the NMSSM stems from either the triple coupling term 
of the Higgs doublets and the Higgs singlet or the cubic term of the Higgs singlet itself, but not 
from the coupling between two Higgs doublets [5-7]. 

In this paper, we impose the restriction of = on the tree-level Higgs potential of the 
NMSSM. We start with the tree-level Higgs potential which is phase-free, by assuming that A, 
A\, k, and Ak are real. Thus, = in our tree-level Higgs potential. Now, the two Higgs 
doublets and the Higgs singlet may be expressed via a unitary transformation as 

vi + Si + i sin pA \ 

cos \ (3) 

V2 + S2 + icos/3A J ' 

N = (x + X + iy), 
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where ^i, X arc the scalar Higgs fields, A, Y are the pseudoscalar Higgs fields, and C+ is 
the charged Higgs field. 

^From the tree-level Higgs potential, the symmetric 3x3 mass matrix Ms for the scalar 
Higgs bosons is obtained in the basis {81,82, X) as 

{Ms)ii = m| cos^ (3 + \x{Ax + kx) tan (3 , 
{Ms)22 = m% sin^ P + >^x{Ax + kx) cot (3 , 

{Ms)33 = {2kxf - kxAk + ^v^^x sill 2/3 , 

{Ms)u = {X^v^ - sill 2/3 - \x{Ax + kx) , 

{Ms)\z = 2\^xvcosP - \vsmf5{Ax + 2kx) , 

{Ms)23 = 2A^OTsin/?- Af cos/3(^A + 2A;x) , (4) 
where the gauge boson masses are given as = [gf + g2)v^ jl and = g^v^ for v = 



' + V2 = 175 GeV after the electroweak symmetry breaking. In the above mass matrix, 
the soft SUSY breaking masses m^j^, "^1^2' ^^"^ already are eliminated by minimization 
conditions with respect to 8\, 82, and X. Thus the mass matrix for the scalar Higgs bosons 
is a function of tan /3 = V2/V1, X, A\, k, Ak and x. The mass matrix M5 is then diagonalized 
in the basis of mass eigenstatcs. Let us denote the mass eigenstates of Ms as Si, 82, 83, which 
are linear combinations of the weak eigenstates, namely, 8i,82,X. (Notations are a little bit 
confusing but hereafter 8i {i = 1,2,3) are the mass eigenstates.) The eigenvalues of the mass 
matrix Ms are obtained and sorted as 

m|. = lTr(M5) + 2VI^cos|^^^| (j = 1, 2, 3) , (5) 



where 



e = cos-M-L) , (6) 



W3 



with 



W = -^{Tr{Ms)}^ + ^TT{MsMs) , 

U = -^{TV(M5)}3 + ^T:t{Ms)Tr{MsMs) + ^ det(M5) . (7) 

Thus, is the squared mass of the mass eigenstate 81, and so on. 

For the pseudoscalar Higgs bosons, the symmetric 2x2 mass matrix Mp is obtained in the 
basis {A,Y) as 



(Mp) 
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2Xx{Ax + kx) 



sin 2/3 

(Mp)22 = 3kxAk + 2Xkv'^ sin 2/3 + — sin 2/3 , 

(Mp)i2 = Xv{Ax - 2kx) . (8) 
The pseudoscalar Higgs boson masses are given by the eigenvalues of Mp as 

^PuP2 = \ [^i^p) T \/(TVMp)2-4det(Mp)| . (9) 
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Finally, the charged Higgs boson mass at the tree level is given as 



2 2 \ 2 2 2Ax . . > 

sin 2(3 



(10) 



Unlike in the case of the MSSM, the tree-level mass of the charged Higgs boson can be either 
heavier or lighter than the W boson according to whether the second is smaller than the third 
term or not. 

Next, we calculate tree-level masses of other relevant fields. These will be used later for 
evaluating the radiative corrections to the neutral Higgs sector. The fcrmion matter fields 
obtain their masses after the clcctrowcak symmetry breaking as = {htV2)'^ for top quark, 

= {hfyVi)^ for bottom quark, and = (hrVi)^ for tau lepton. The scalar top quarks, the 
scalar bottom quarks, the scalar tau leptons, and the charginos have their tree-level masses as 



tl,2 



mt + ^irnl + m\) + ^ cos 2/3 ^ 



+ 1 \mw 



— ) cos 2/? 



01,2 



+ mf {Af + X^x^ cot^ /? + 2AtXx cot /? cos (t)t) 

ml + l {ml + m|) - ^ cos 2/? ^ 



+ 



1 



m| 



1 



-m 



W 



12 ^ 3 



COS 2p 



+ mi{Af + X^x'^ tan^ /3 + 2Ai,Xx tan (3 cos ( 



ml 



ml + \{ml + ml) - ^ cos 2(3 ^ 



|^(mi -m|) 



1 

-m 



w 



cos 2(3 



+ mi{A^ + X^x"^ tan^ (3 + 2ArXx tan /3 cos ( 



1 



Xl,2 



I ^ (M| - A V) - cos 2/3 y 



-I- 2m%^ cos^ (3{Ml >?x^ tan^ /3 -|- 2M2 Aa; tan (3 cos ^c) 



(11) 



where mg, my, ms, mi, and mg are the soft SUSY breaking breaking masses, and A^^ Af,, and 
At are the trilinear soft SUSY breaking parameters of mass dimension. One can find that there 
are four phases (f>t, (pbi (pr, and 0c in the above expressions: They are originated from the trilinear 
soft SUSY breaking parameters At, A),, At and the SU(2) gaugino mass M2 which are complex 
in general whereas A is assumed to be real and positive. Note that the scalar fermion masses are 
not symmetric under mq my, mq ^ ms, or mi, <r^ mE since D-terms are included in them. 
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In the NMSSM, there are five neutrahnos. Their tree-level masses are obtained from the 5 
X 5 mass matrix 












\ 





—Mo 




V2 ' 





V2 ' 







XN 




-32- H'^ 
V2 ' 


92 ffO 


XN 





XHl 












-2kN ) 



(12) 



V 

where Mi is the U(l) gaugino mass. Note that we have already assumed that A and k are real. 
Thus, the ncutralino mass matrix may have two complex numbers, M\ and Mi. Let the relative 
phase between them be This is therefore the only non-trivial phase in the neutralino mass 
matrix in our case. 

The above neutralino mass matrix is complex and symmetric, but not Hermitian. In other 



to calculate the neutralino masses, the Hermitian matrix A^^oAlto is diagonalized through a 

,2 



similarity transformation. The squared neutralino masses at the tree level are denoted as m^o 
(i = 1 to 5), and they are sorted such that m-o < m-q for % < j. If they are all different from 



each other, they are obtained as the solutions of 

which may be rewritten as 

ml^, + im|o + Bmlo + Cmt, + Dmlo + E 

Xi Xi Xi Xi Xi 







(13) 
(14) 



where the field-dependent coefficients A, B, C, D, and E are very complicated and obtainable 
from Appendix A. 



III. Radiative corrections 

Now let us consider the radiative corrections. In order to describe the explicit CP violation 
scenario, it would be convenient to decompose the neutral Higgs fields as 

Hi = ( vi + hi + i sin ph^ ^ , 

= (v2 + h2 + icosf3h3^, (15) 
N = y X + h4 + ih^ j , 

where hi {i = 1 to 5) are the neutral Higgs fields which do not have definite CP parities. In our 
case with explicit CP violation scenario, they are mixtures of CP-even and CP-odd components. 
Through radiative corrections, there emerges CP mixing among hi. On the other hand, the 
charged Higgs bosons does not mix because both of them are CP-even. 

Now, we turn to the radiative corrections. The incomplete cancellation between ordinary 
particles and their superpartners give the one-loop corrections to the tree-level Higgs boson 



6 



masses. The full Higgs potential up to the one-loop level is the sum of the tree-level Higgs 
potential and the radiative corrections as 

where may be decomposed into 

= 1/* + yf' + + + yx° ^ (16) 

with superscripts indicate the radiative corrections due to the relevant particles and their super- 
partners. In fact, the neutral part and the charged one are related to each other. In the NMSSM 
with explicit CP violation, all the neutral Higgs boson masses can be expressed in terms of the 
charged Higgs boson mass. Note that the charged Higgs boson mass should be used as input. 
And and should be evaluated simultaneously. 

We employ the effective potential method in order to calculate radiative corrections [10]. 
The general expression for the one-loop effective potential is given by 

= E ^(-l)''H2J. + l){M,ihO}' 

k 

where Jk is the spin of particle or superparticle, A is the renormalization scale, and Cfc = 
CcoiorCcharge with Ccoior being the color factor and Ccharge the charge factor of the corresponding 
particle or superparticle. We have Ccoior = 3 (1) for colored (uncolored) particles and Ccharge = 
2 (1) for charged (neutral) particles. Note that depends on the neutral Higgs fields hi {i = 
1 to 5). 

Explicitly, we write down the formulae for the relevant one-loop effective potentials: The 
one-loop effective potential for the contributions of top quark and scalar top quarks is given by 

the one for bottom and scalar bottom quarks by 

the one tau lepton and scalar tau leptons by 

^ 32i2 A2 2) 16n^ \ ^ 2) • 

the one for W boson, the charged Higgs boson, and the charginos by 

= ^ ^ - 2 J + 3^ ^ - 2 J - g ^ - 2 J ' 

and lastly, the one for Z boson, the neutral Higgs bosons and the neutralinos is given by 
_.,c 2.M%(. Ml 3\'Ml( Ml 3\ ^^h(, A 



log 



{Mkjhj)}' 
A2 
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Note that in the above one-loop effective potentials fermions enter with a negative sign while 
bosons with a positive sign. 

If CP be conserved in the Higgs sector, /13 and /15 would be the pseudoscalar components of 
the Higgs fields, and the minimum conditions for them would be trivial. In our case of explicit 
CP violation, the minimum conditions for them are non-trivial. Thus, we obtain two CP-odd 
tadpole minimum conditions for and /15, respectively, as 

Sm^AtXx sin (j)t 2 2^ , 3m? AAx sin 2 2x 

Wn^v^sm^P ^ *i 167r2w2 cos2 /? &i' 

m?-AT-Xx siiKpT .,0 2 \ myi^XxM2sm(f)c „/ 2 2 % 
-/(m^^,m^) ^-^^ /Ki'^xa) 




Bzmlo + Csm^o + Dsmlo + E3 
n (m|o -m|o) 



(17) 



3m(2^tAsin(/>t 2 2^ , Sm^AAsin^;, 2 2n 

= — -T. o — f{m~,mf]-\ ^ TTTT fi'iTT'i ,mi ) 

167r2t;sin2/3 ■'^ *i' 167r2t; cos2 /? ^ ^-i' 

m2^rAsin0^ 2 2n , rnfj^AM2 sin (/^^ 2 2 x 

+ le^^w^a -^^"^-'^-^ + — — ^Ki'^xJ 

+ ^5il67r2^;sin2/3|^°^ j " ^ j ^^K^^^^?^ ' ^^^^ 

where the five terms come respectively from the scalar top quark, the scalar bottom quark, the 
scalar tau lepton, the chargino, and the neutralino contributions. The dimensionless function 
arising from radiative corrections is defined by 

1 1 7?7'2 TnP 

and we need to introduce Aj = {dA/dhi), Bi = {dB/dhi), Ci = (dC/dhi), Di = (dD/dhi), and 
£'i = ipE/dhi): In the above minimum conditions, just and E'j (i = 3 and 5) enter. 

Note that no term associated with the tree-level is present in the above minimum conditions 
since (j) = Q. 

By differentiating the Higgs potential V at the one-loop level with respect to the neutral 
Higgs fields, a 5 x 5 symmetric mass matrix M for them is obtained in the {hi,h2,hs,h4,h5)- 
basis, which may be decomposed as 

where M-^- is obtained from while Mlj from V^. Explicitly, M^j may be expressed in terms 
of the scalar and pseudoscalar mass matrix elements as 






= iMs: 


hi 


M2°2 


= {Ms: 


)22 


M44 


= {Ms: 


)33 


M?2 


= {Ms, 


)12 


M?4 


= {Ms, 


)13 
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= {Ms)23 

and 

= (Mp)ii 

M°5 = (Mp)22 
M°5 = (Mp)i2 

and the remaining elements of M^j are explicitly zero since we assume that = 0. 

Thus, the elements of the neutral Higgs boson mass matrix at the one-loop level may be 
written as 

Mil = (M5) 11 + m\ sin^ (3 + Mi\ , 

M22 = (M5)22 + ml cos^ (3 + M2^2 , 

M33 = (Mp)ii +m\ + M3I3 , 

M44 = (M5)33 + m| + Afii , 

M55 = (Mp)22 + m| + , 

M12 = (M5) 12 - m\ cos /3 sin /3 + Mi^2 > 

Ml 3 = M^ 



'13 = -(Wis ) 
Mi4 = (M5)l3-^^+Mi 



-y cos/3 



14 ! 



Mi5 = MI5, 
M23 = M2^3 , 

M24 = (M5)23-^^+Ml4, 

M25 = M2I5 , 
M34 = M3I4 , 

N rriiv COS 3 sin 6 ,1 

M35 = (Mp)i2 + ^ + M3I5 , 

M45 = Mis , (19) 
where two expressions are introduced for convenience as 



2 _ SnifAtXx cos (pt 2 2 X 3m^At,Aa:cos(/)b 2 2n 
~ 167r2t;2 sin3 /3 cos /3 ^ ^""^i ' ""^^ ' IQn^v^ cos^ /3 sin /3 ^ ^""^i ' ' 



ml AtXx cos 4>r 2 2 x m^\xM2 COS (t)c 2 2x 

^ 167r2^;2 cos3 psinp^ ^""^^ ' ""^ ^ 2Tr^v'' sin 2/3 ^""^i ' ""^^ ^ 

5 r / \ ] i33m^o + ^33m?o + Cssmio + .033^20 + E33 



and 



^ 3mf^tAcos0t 2 2n , 3mgAbAcos0b 2 2^ 

'TIr = 7i 7^t{mf,mf)-\ ^ 7^f{mt ,mi , 

^ 167r2xtan/3 ^ *i' 167r2xcot/3 ^ ''1 

,2 /I \^^^A ^2 



m^ArAcos^T,. 2 2\ m4.AM2sin2/3cos0c 2 2 x 
+ 167r2.cot/3 ^("^-'"^-) §^^2^ /K.'-x.) 

5 rn|o f / mio \ ] A^^mlo + B^bm%o + C^^mt, + D^^mio + ^55 

I \ Afc ) 1 „ ''^fc Ilk ±k ±k Afe Xfc 
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In each expression, the first term comes from the scalar top quark contributions, the second 
term from the scalar bottom quarks, the third term from the scalar tau Icptons, the fourth one 
from charginos, and the last term from neutralinos. In niA and mg, the coefficients A^, Ba, 
Cii, Da, and Eii (i = 1 to 5) are given by Aij = \/ijA, Bij = SJijB, = ^ijC, Dij = SJijD, 
and Eij = \/ijE, where the second-order differential operator is defined as 

d'^ ^ d ^ _]__d_g / sin^/3 d cos^ ^ \ s 

dhidhj vi dhi V2 dh2 \ vi dhi V2 dh2 J 
1 ^ X ^ 

X dh4 X dh4 ' ^'^^^ 

for j >% = 1 to 5). In the above differential operator, all the terms except the first one are 
caused by the shift in the vacuum due to the radiative corrections. 
We may further decompose Mlj as 

where M-j is obtained from F*, M^j from V^, MJ^ from , Mf- from V^, and M|° from , 
after imposing the CP-odd tadpole minimum conditions for ^3 and h^. We note that Mf- and 
Mf- are matrices for squared masses. The elements of M'^ have no color factor. In Appendix 
B are presented the complicated formulae for Mly In Appendix C for M^^-, in Appendix D for 
MJ^, and in Appendix E for Mf-. 

Now, we concentrate on M^- , which represents the radiative corrections due to the contribu- 
tions from the neutralino sector: It consists of the contributions from Z boson, tree-level scalar 
and pseudoscalar Higgs bosons, and neutralinos. Thus, we may decompose it further again as 

= A^. + A^. + A^ + A^° 

where each term represents the relevant contribution. Let us describe one by one. 

^0 ^0 " -0 -0 
The contribution of Z boson affects the four elements Mj\ , M22 , and Mj^2 = -^21 only, 

since the Higgs singlet A'' does not couple to Z boson. Calculations yield the contribution of Z 

boson as 

aZ \ / „„^2 o \ 

^log ( , (21) 





/ COs2 p \ 




sin2 p 




\^ cos/? sin/? J 



16^ ^Va 



with Afi = Afo. 



*2i — ■"12 ■ 

The contribution of the tree-level scalar Higgs bosons is obtained as 




*j 1=1 

where {I = 1, 2, 3) are the tree-level masses of the scalar Higgs bosons of Eq. (5), and 

,2 A 4 , D_2 



dhi {rnl^ - - mlj ' 
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and 



dhidhj (m|^ - m\ ) (m| ^ - ) 



1 / dm% dm% dm% dmi \ , , 

(m|^ — m|^) I 50j d^i d4>j J ' 

for i = 1 to 5, and / rrig^ 7^ "^5c' -'^^'^^j ^® introduce 

A = - Tt{Ms) , 

B = Ms^^Ms22 + Ms,,Msss + MS22MSSS - Ml^, - - mI^^ , 

C = - det{Ms) , (25) 

where Msij are given in Appendix F, and Ai = (dA/dhi), Bi = (dB/dhi), and Ci = {dC/dhi). 

As aforementioned, all terms in Vij except the first one are caused by the shifts in the 
vacuum due to the radiative corrections. Note that the shifts take place along /ii, /i2; and /14. 
Those shifts occur actually on only the diagonal elements of the mass matrix of the neutral 
Higgs boson through Vij: Mn by the second term of V) -^22 by the third term, M33 by the 
fourth term, M44 by the fifth term, and M55 by the sixth term, respectively. Since the tree-level 
Higgs potential does not possess the scalar-pseudoscalar mixing, both and become the 
pscudoscalar Higgs fields and dV^/dh^ = dV^/dh^, = 0. At the one-loop level the non-trivial 
solutions for and already are obtained by the CP-odd tadpole minimum conditions. 

Next, the contribution of the tree-level pseudoscalar Higgs bosons is obtained as 

1 ( dApdAp dApdWp \ log{mlJml^) {\/ijWp) 2 2^ 

647r2 [dhi dhj ^ dhj dhi ) {ml^-mlj Uti^ ^^""^i' ""^^i 

1 ( dWpdWp \ g{m\, ml;) 1 ( dApdAp \ ( mlml \ 

647r2 1^ dhi dhj ) {m\ - mj,J^ 6An^ \ dhi dhj ) \ A^ J ^ ^ 

for i,j = l to 5, and mp^ and mp^ are the tree-level masses of pseudoscalar Higgs bosons. Here, 
we have introduced Ap and Wp as 




A 



p 



iMp,,+Mp,,) 



2 

\2 



Wp = (^^--^^-^ +Ml,^. (27) 



where -Mpi^ are given in Appendix G. 

Lastly, the contribution of the neutralinos is obtained as 

AS=-g^ logj^-l U^-g5^1o.^ U^ . (28) 



where 

dml Aivnl + Bivnl + Qm^ + Diml + % 



dhi {ml -ml) {ml - ml ) {ml - ml) {ml - ml ) 



(29) 
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and 



dhidhj (m^ — m^)(m^ — m^)(m^ — m^)(m^ — m^) 

+ 



1 / dm^ dml dml dm'^ 



{ml — ml) \ dhi dhj dhi dhj I 

1 / dm'^ dml ^"^c ^"^a \ 
(m^ — m^) I dhj J 

^ 1 dm^dm^ ^ dm^dm^\ 

{ml — m§) y dhi dhj dhi dhj J 

^ 1 / dml dml ^ dml dml \ 
{ml — ml) \ dhi dhj dhi dhj J 

for i,j = l to 5. In the above two expressions, we use the shorthand notations for the neutrahno 
masses: m„ = m|o and so on. We assume that none of the neutrahno masses are equal one 
another. 

Wc note that Afj in above equation is not in its final form, in the sense that the two CP- 
odd tadpole minimum conditions should be applied to obtain the explicit results. The minimum 
conditions on hi, h2, and /i4 are already imposed through the differential operators Vij. 



IV. Numerical Results 

Experimental lower bounds on the Higgs boson masses in various models have been set by the 
LEP2 data. They have excluded at the 95 % confidence level a Higgs boson in the Standard 
Model, which is lighter than 114.4 GeV [11]. For the MSSM without explicit CP violation, the 
lower bounds on the scalar and pscudoscalar Higgs boson masses have been set as about 84 GeV 
and 86 GeV, respectively, by the data [12]. For the MSSM with explicit CP violation, on the 
other hand, the data have allowed the possibility of the lightest neutral Higgs boson with a mass 
of 60-70 GeV [13]. 

Compared to the MSSM, the Higgs sector of the NMSSM have seemed not utilized the LEP2 
data in full. Some years ago, the analysis of the LEP2 data have shown that a masslcss scalar 
Higgs boson might exist in the NMSSM without explicit CP violation [9]. This analysis have 
been done when the center of mass energy and the luminosity of the LEP2 have not reached their 
final capacity [9]. In the NMSSM with explicit CP violation, the mass of the lightest neutral 
Higgs boson has been roughly evaluated with LEP2 data [6]. 

If CP symmetry be conserved in the Higgs sector of the NMSSM, the five Higgs fields 
decouple such that hi, /i2, and /14 would become the scalar Higgs bosons while hs and ^5 the 
pseudoscalar Higgs bosons. There would be no mixing between them. The CP violation is 
dictated by the six off-diagonal elements of the mass matrix, M13, M23, M34, M15, M25, and 
M45, of the neutral Higgs bosons in the basis {hi,h2,h^,hii,h^), which describe the scalar- 
pseudoscalar mixing. The remaining off-diagonal elements describe the scalar-scalar mixing or 
the pseudoscalar-pseudoscalar mixing, and have nothing to do with the CP violation. 

For numerical analysis, we set the fermion masses as mt = 175 GeV, mj, = 4 GeV, and m-T = 
1.7 GeV, the gauge boson masses as mz = 91.1 GcV and mw = 80.4 GeV, and the weak mixing 
angle as sin^ Ow = 0.23. By assuming universal gaugino masses at GUT scale, the low-energy 
values of the U(l) and SU(2) gaugino masses has a mass relation of Mi = 5/3 tan^ 9w^2- We 
take the assumption. 
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Now, to be specific, we take a representative point in the parameter space: We choose tan (3 
= 5, A = 0.5, k = 0.4, Ax = 590 GeV, Ak = 5 GeV, x = 178 GcV, ttiq = rriL = 967 GcV, 
ruT = ruB = THE = 79A GeV, At = At, = A^- = 907 GeV, and M2 = 357 GeV. The complex 
phases except the one from neutralino contributions are fixed as (pt = 4>b = 4>t = = IOtt/IQ. 
Then, we are left essentially with two free parameters: A and ^1. The phase ^1 from the 
neutralino contributions is allowed to vary as < </>i < 27r. 

We would like to mention about the dependence on the renormalization scale of the neutral 
Higgs boson masses. Physically observable quantities such as the neutral Higgs boson masses 
should be scale-independent. However, Fig. 1 shows a weak scale-dependence of the neutral 
Higgs boson masses, where m,h^ [i = 1 to 5) are plotted against A for (pi = 107r/19, and the 
other parameter values are set as above. This is because we consider only up to the one-loop 
correction: If all of the higher-order contributions are further taken into account, the scale- 
dependence would disappear. Hereafter, in our analysis, we fix A as 300 GeV in the one-loop 
effective potential. 

We are interested in the effect of on the 5x5 neutral Higgs boson mass matrix M and 
the 5x5 orthogonal mass matrix O that diagonalizes it. Among the 10 off-diagonal elements 
of the symmetric M, six of them describe the scalar-pseudoscalar mixing: M13, M23, -M34, M15, 
M25, and M45. The scalar-pseudoscalar mixing, which is responsible for CP violation in the 
NMSSM, is triggered by five phases in these matrix elements, namely, ^j, 0,-, and (f)i. If 
sin((^t) = sin((?!)b) = sin(0T-) = sin((^c) = sin(0i) = 0, these matrix elements would be zero, and 
the scalar-pseudoscalar mixings are said to be minimal. Thus, the scalar-pseudoscalar mixings 
would disappear to restore the CP symmetry, if the Higgs sector contains no complex phases. 
In this case, O would also be simplified as O13, O23, O34, O15, 025, and O45 would be zero. 

In Fig. 2, we first plot the numerical result for the neutral Higgs boson masses, as functions 
of 01, where A = 300 GeV and the values of the other parameters are the same as Fig. 1. The 
mass of the lightest neutral Higgs boson is about 80 GeV, and roughly stable, for our choice of 
parameter values. The mass of the next-to-lightest neutral Higgs boson, on the other hand, is 
found to vary by 40 % for 01 between and 27r. This indicates that the neutralino contributions 
are not negligible for some parameter values. 

We then plot the dependence on 0i of M13, M23, M34, M15, M25, and M45 in Figs. 3a and 
3b, for the same parameter values as Fig. 2. The unit is (GeV)^. In Fig. 3a, three curves 
for Mi3 (solid curve), M23 (dashed curve), and M34 (dotted curve) are displayed, while in Fig. 
3b those for M15 (solid curve), M25 (dashed curve), M45 (dotted curve). It can be seen that 
the scalar-pseudoscalar mixing elements for the neutral Higgs bosons fluctuate widely in the 
NMSSM with explicit CP violation, by the one-loop neutralino contributions. 

These results may be compared to the case of the MSSM. In the case of the MSSM with 
explicit CP violation, there is a complex phase and only two scalar-pseudoscalar mixing elements. 
As the relevant complex phase varies between and 27r, it is found that one of the scalar- 
pseudoscalar mixing elements in the MSSM increases as much as twice its value at zero complex 
phase, due to the neutralino contribution [3]. 

In order to obtain a more comprehensive picture on the amount of mixing between scalar 
Higgs bosons and pseudoscalar Higgs bosons, which is directly linked to the amount of CP 
violation the the Higgs sector of the NMSSM, we would not only examine the magnitudes of 
Mi3, M23, M34, Mi5, M25, and M45, but also examine the relative sizes of them against the other 
matrix elements. The ratios of the mixing matrix elements to the non-mixing matrix elements 
would tell us the relative importance of the scalar-pseudoscalar mixing elements. 
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We define the ratios of mixing matrix elements to non-mixing ones by 



where are (1,3), (2,3), (3,4), (1,5), (2,5), and (4,5) and {p,q) are (1,2), (1,4), (2,4), and 
(3,5). The number of matrix elements that do not mix the CP states is 4: M12, M14, M24, and 
M35. Since we have 6 mixing matrix elements, the number of possible ratios in the NMSSM is 
24, whereas in the MSSM there are just two ratios. 

In Figs. 4a to 4h, we plot all of these 24 ratios of the mixing matrix elements to the non- 
mixing ones, as functions of (pi, for the same parameter values as Fig. 2. We check that the 
absolute values of non-mixing elements M12, M14, M24, and M35 are much larger than 1 (GeV)^ 
for Figs. 4. These ratios vary largely as 0i varies from to 27r. On all Figs. 4a to 4h, these 
ratios reach their peak values at the values of 01 where the mixing matrix elements themselves 
become extremum values. ^From Figs. 4, one may deduce that for a particular set of parameter 
values the scalar-pseudoscalar mixing matrix elements are comparable in size to the non-mixing 
matrix elements. 

On the other hand, we can randomly explore the parameter space to check the values of 

(i i) 

ip<i)' '^^^^ because the results in Figs. 4 might not faithfully represent the whole parameter 
space of the NMSSM, since, up to now, our analysis is carried out for a specific set of parameter 
values (except ^i). Thus, we would like to search in more detail the meaningful region of the 
parameter space in order to see if the mass matrix elements can be still comparable in size there. 
The boundaries of the region of the parameter space we search are defined by 

A = 300 GeV 

< (f)t = (f>b = (f>r = (l>c < '^'TT 

2 < tan /? < 40 
< A, A: < 0.7 

< Ax,Ak,x,M2 < 1000 GeV 
0<At=Ab = Ar< 1000 GeV 
0<mQ = mL< 1000 GeV 
0<mT = mB = mE< 1000 GeV 



and 



< 01 < 27r . 



(i i) 

We evaluate g) at a point in this region, which represent naturally a given set of parameter 

values. We take three values of 0i: 01 = 0, 0i = 7r/3, and 0i = 27r/3. For each value of 0i, we 

(i i) 

choose 20,000 points by Monte Carlo method in the above region, and at each point 24 

(i i) 

are evaluated. Thus, for each value of 0i, we obtain 480,000 numbers for R\^ qy 

We put these 480,000 numbers into four categories for convenience: -R(^'"'g') > 1, 1 > -^(pg") ^ 

0.1, 0.1 > R^l^l^ > 0.01, and 0.01 > R'j^^^y In Table 1, we count them and list the result. One 

(i i) 

can see in Table 1 that, for example, for 0i = 0, 188 of 480,000 R(^ q) are larger than 1. 

(i i) 

It can roughly be said that the larger g) are the larger mixings between CP states become. 
^From Table 1, it is apparent that we have clearly more large -^1*' for 0i 7^ than for 0i = 0. 
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Table 1: For each of (pi = 0, vr/3, and 27r/3, we evaluate 24 at 20,000 points in the region 

of the parameter space defined in the text. The 480,000 values of -^(p are classified into four 
categories. We list the counts in each category. 





B}i'''\ > 1 
ip,q) - 


1 > <'J) ^ 


0.1 > Rf;% > 0.01 


0-01 > 


</.i = 


188 


1,525 


13,614 


464,673 


(pi = 7r/3 


1,051 


9,056 


66,312 


403,581 


(pi = 27r/3 


1,417 


12,807 


73,584 


392,192 



The number of large -R(pg^)5 larger than 0.1, increases as (pi changes from to 27r/3. This may 
be interpreted as the effect of (pi, the complex phase in the neutralino sector. 

Now, we turn our attention to the 5x5 orthogonal matrix that diagonalizes the mass matrix 
of the neutral Higgs bosons. This orthogonal matrix also contains pieces of information about 
the CP violation and the mixing between scalar and pseudoscalar Higgs bosons. The elements 
Oij determine the couplings of the physical neutral Higgs bosons to other fields in the NMSSM. 
In Figs. 5a and 5b, we plot as a function of (pi, for the same parameter values as Fig. 2, the 
squares of O13, O23, O34, O15, O25, and O45, which link the CP-even and CP-odd states. One can 
see that O34, Of 5, and O25 fluctuate widely for < ^1 < 27r, while the remaining ones remain 
nearly zero. Therefore, for the particular set of parameter values, it is possible to expect that 
neutralino contributions may affect seriously the scalar-pseudoscalar mixings in the NMSSM. 
Note that, in our analysis, ^1 = does not imply the absence of CP violation, because there are 
other non-zero complex phases, i.e., (pt = (pb = (pr = (pc = 107r/19. 

Meanwhile, one can evaluate effectively the size of the CP violation by measuring the di- 
mensionless parameter [6] 

P = 5'OnOiiOiiO|iOii , 

The range of p is from to 1 since the elements of the transformation matrix satisfy the orthog- 
onality condition of X]j=i Oji = 1. If p = 0, there is no CP violation in the Higgs sector of the 
NMSSM. On the other hand, if p = 1, CP symmetry is maximally violated. The maximal CP 
violation that leads to p = 1 takes place when Of^ = Ofi = O31 = O41 = Of^ = 1/5. In Fig. 6, 
we plot p, as a function of (pi, for the same parameter values as Fig. 2. One can observe that 
p reaches to ~ 10~^ for certain values of (pi. This shows that it is possible for the neutralino 
sector to play a dominant role in the CP violation in the NMSSM at the one-loop level, through 
the scalar-pseudoscalar mixings, at least for the particular set of parameter values. 

Here, too, we explore the parameter space randomly and check the value of p. The regions 

(i i) 

of the parameter space is the same as the aforementioned one for R(^p qy We take 19 values of 
(pi from to 27r. For each value of (pi, we choose 20,000 points in the region and evaluate p, and 
take the largest value among them. The result is shown in Fig. 7, where the maximum value 
of p, say praax, for each (pi is plotted. We see that pmax is considerably large for a wide range 
of (pi. Thus, the CP violation can be significantly affected by the presence of (pi. Without the 
neutralino contributions in the NMSSM, the explicit CP violation at the one-loop level might 
be considerably reduced. 

We now estimate the dependence on (pi of the mass of the lightest neutral Higgs boson m/^ 
in more detail. We first choose randomly a point in the region of the parameter space defined 
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above. We check at this point if the condition of — i^h(^„) > 10 GeV, (n = 1 to 4) is 

satisfied. Then, we change only the value of (pi arbitrarily while fixing other parameter values 
of the point: This defines another point. We check at the new point too if the condition of 
™'i(n+i) ~ "^/^cn) > lO GeV, (n = 1 to 4) is satisfied. If so, we keep the two points as a pair. Let 
the values of (pi of a given pair of the two points be ci and C2. 

Now, at each of the the two points, we evaluate the mass of the lightest neutral Higgs boson, 
m/ij , and introduce the ratio of the mass values as 

^ ^ mhi(.(t>i = ci) 
mhi{(pi = C2) 

for mhi{(pi = ci) > mh^{(j)i = C2), and 

^ mh^{4'i = C2) 
mhi{(pi = ci) 

for mhi{(pi = ci) < mh^{(pi = C2). In this way, we always have Ri > 1. Note that Ri = 1 
implies the independence of the lightest neutral Higgs boson mass on (pi, i.e., on the neutralino 
contributions. 

We select 20, 000 pairs of points randomly in the aforementioned region of the parameter 
space and evaluate Ri for each pair. We plot the 20,000 points of Ri against C3 = |c2 — ci| in 
Fig. 8. One can see in Fig. 8 that only a few of the points are scattered away from the line of 
Ri = 1. This indicates that for most of the parameter region the (pi dependence of the lightest 
neutral Higgs boson mass is weak. In other words, the neutralino contributions to the lightest 
neutral Higgs boson mass in the NMSSM are rather small. The masses of other neutral Higgs 
bosons behave similarly. Thus, one might say that the masses of the five neutral Higgs bosons 
are roughly independent on the CP phase of the neutralino sector. 



V. Conclusions 

We have investigated the possibility of explicit CP violation at the one-loop level in the neutral 
Higgs sector of the NMSSM, where the tree-level Higgs potential is assumed to possess CP 
symmetry. The mixing between the scalar and pscudoscalar Higgs bosons, which induce the CP 
violation, occurs via the radiative corrections at the one-loop level. For the radiative corrections, 
we consider the loop contributions due to top quark, bottom quark, tau lepton, W boson, Z 
boson, the charged Higgs boson, five the neutral Higgs bosons and their superpartners. With 
those contributions, five phases appear. They are: (pt from the scalar top quark sector, (ph from 
the scalar bottom quark sector, (pT from the tau lepton sector, (pc from the chargino sector, and 
01 from the neutralino sector. 

The effects of CP violation in the NMSSM at the one-loop level may manifest themselves 
in many places: The mass matrix of the neutral Higgs bosons contain some elements that mix 
CP-even and CP-odd fields. Some elements of the orthogonal matrix that diagonalizes the mass 
matrix of the neutral Higgs bosons have also carry information about the CP violation. They 
will be zero if CP is conserved. The five complex phases are all responsible for, and contribute 
to, the CP violation. We pay attention to (pi, and examine the size of contributions from the 
neutralino sector to the matrix elements relevant for the CP violation in the NMSSM at the 
one-loop level. 

As (pi is allowed to vary from to 27r while the other complex phases are fixed, we can see 
the dependence on <p)i of the CP violation. Prom Figs. 2 to Fig. 6, we observe that for a certain 
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set of parameter values most of the relevant quantities that we examine depend significantly on 
One of the scalar-pseudoscalar mixing matrix elements change from nearly zero to as large 
as more than 10,000 (GeV)^, depending on cpi. Actually, the absolute values of these scalar- 
pseudoscalar mixing matrix elements do not convey more useful information than the relative 
sizes among the mixing matrix elements and non-mixing matrix elements. 

The ratios of the mixing matrix elements and non-mixing matrix elements are examined. 
They also depend clearly on (pi, exhibiting the neutralino contributions. Some of them approach 
to 1, for certain value of 01, and many of them are larger than 0.01 for a wide range of The 
results imply that the scalar-pseudoscalar mixing matrix elements may be comparable in sizes 
to the non-mixing matrix elements, for a certain set of parameter values in the NMSSM, due 
to the neutralino contributions at the one-loop level. The masses themselves of the five neutral 
Higgs bosons are not so affected by 4>i, at least for the certain set of parameter values. 

The elements of the orthogonal matrix that diagonalizes the mass matrix of the neutral 
Higgs bosons carry more direct information on the CP violation than the mass matrix itself. As 
can be seen in Figs. 5, some of the matrix elements Ofj fluctuate very widely between and 1, 
depending on 0i. Moreover, p, which is a quantity that measures directly the amount of the CP 
violation, is also checked. It is found that p is larger than 0.01 for some values of 0i, whereas 
p is nearly zero when 0i = 0. Thus, it is reasonable to expect that the neutralino sector at 
the one-loop level dominantly contributes the CP violation in the NMSSM, for a certain set of 
parameter values. 

We also explore a meaningful region of the parameter space by Monte Carlo method in order 
to check the above observations are not confined to the particular set of parameter values. In 
the parameter region, we pick randomly 20,000 points for given (pi to evaluate p at them. We 
choose the largest of them and denote it pmax for given (pi. It is found that pmax shows a clear 
dependence on (pi. Moreover, pmax at 01 = is definitely smaller than 

Pmax at non-zero (pi 

values. Actually, for some values of (pi, Pmax close to as large as 0.29. To conclude, we find that 
for some parameter values in the NMSSM with explicit CP violation at the one-loop level the 
neutralino contributions to the scalar-pseudoscalar mixing exhibit relatively significant effects. 
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Appendix A 



The elements for the field-dependent mass matrix M.yoM.\o of the neutralinos at the tree level 

^ X, 

are 





= 
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2 

|(|ifOp + |i/Op)+M|, 




= 






= 






= 




A12 
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(l^iT + l^2T), 


X13 


= 




A 14 


= 




A15 


= 




A23 
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--^ (iff M2 + AAr*i^^) , 


A 24 




4 (i?r M2 + AiV*iJ?) , 


A25 




^A(i/?<-<i/o), 


A34 






A35 




A (AiVi^r - 2kN*Hl) , 


A45 




A (AiVi/f - 2kN*Hfj . 


Appendix B 





(31) 



The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the top quark and scalar top quarks are 
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Appendix C 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the bottom quark and scalar bottom quarks are 
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Appendix D 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the tau lepton and scalar tau leptons are 

Ml, = 




22 



sin^ P 



+ 



m 



I sin /3 



4 

rniXxA 



T2 



sin/3A^l log (m? /m? 



cos/? 



m^A^x^^^ sin^ 0r sI^t- 



Stt^u^ cos^ /? (mlj — m: 
Stt^ cot^ (3 {rrif — rUf )^ 



T2 ""ny 



Stt^ cot^ /3 (m: 



m^AxAfj 
cos /3 



- 

sin /3Af ' 



cos P 



+ 



cos/jA^ I g(m|^, m|) 
2 /(m?-m?)2 



+ 



+ 



'3m| 
4 ' 

Ami 



+ 



sin 2/3 
32^ 

cos /3 

327r2v2 \^cos2/3 
m| sin 2/3 j ni^AT-Af^ 



+ 



cos^ /3 

4777.2 



647r27;2 

777^ sin 2/3 
2567r27;2 I cos2/3 

7772Ax^t- sin (j!), 

STr^T;^ cos2 /3 
m^XxAT sin 




7772AxAf2 

cos/3 



sin/?Af 1 \og{mf^/mf 



{m: 



,) 



10g("^f2/"^fl) 



777? J 



7n,2 A^-Af-^^ 



Wtt'^v^ cos /3 I cos2 /3 



m'^XAf^ 
Stt'^v cot /3 



rn^ At Af^ ^ cos/3Af 




cos /3 

2777,2 



i?(r77? 777? 



+ 




7772 A COS /3Af, 
327r27; cot /3 

vo^XAt sin 0^ 
87r27; cot /3 

rn^XAr cos /3 sin 
327r27;cot/3 

7772 Ax^T- sin / 7772AxAf2 sin/3Af 
87r27;2 cos2 /3 \^ cos^ 2 

m^rri^XxAT tan /3 sin log(r772^ / m~._^ ) 
327r27;2 cos /3 (777?^ — tti?^ ) 

sin/3A 



777^^ - 777t^, 
log(777? /777? 



2 _ 
^2 

5("7|i; 



777J 



7/7, 



2 ^2 



mlXAf^ 



m^XxAf^ 



87r2'ucot/3 \ cos/3 



fJ 7772 A2x tan /? 
(777? — 777? )2 87r27; cos /3 



/(777: 



-Tl' 



23 



M34 



where 



+ 



327r2?;cot/3 



(m? — 



r2 Ti^ 

25 ~ 87r2t;cot/3 cos/3 2 ) (m\-m\f 

mlmzArX sin /3 sin (pr ^og{m'^^/mU 
32'k'^v cot 6 (m?- — rn^ ) ' 

^ m^A^xA^sin^^Afj ^("^li, "^Ij 



Stt^i; cos^ /3 cot /3 (ml^ — )2 



^ _ mfX^xA^ sin^ ff(m|^ , m|J 
^5 ~ STT^i; cos2 /3 cot p (m?^ - J' 



^5 87r2cot2/3 (rn2^_rn|J2' l^^j 



Af^ = At + Xx tan /3 cos , 

Afj = cos (f>r + A.T tan /? , 

= l^-m^l |mi-m|+|^-m^)cos2/3^ . (37) 




Appendix E 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the W boson, charged Higgs boson, and charginos are 
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Appendix F 

The elements for the field-dependent mass matrix of the scalar Higgs bosons at the tree level 
are 
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The elements for the field-dependent mass matrix of the pseudoscalar Higgs bosons at the tree 
level are 
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Figure Captions 



Fig. 1 : The plot of the neutral Higgs boson masses as functions of A for (pi = IOtt/IQ. We 
set tan/3 = 5, A = 0.5, k = 0.4, Ax = 590 GeV, At = 5 GeV, x = 178 GeV, niQ = mi 
= 967 GeV, mr = niB = tue = 794 GeV, At = A^ = A^ = 907 GeV, Ma = 357 GeV, 
4>t = ^h = = ^c = 107r/19. 

Fig. 2 : The plot of the neutral Higgs boson masses as functions of 01 for A = 300 GeV. The 
remaining parameters have the same values as in Fig. 1. 

Fig. 3a : The plot of mass matrix elements M13 (solid curve), M23 (dashed curve), and M34 
(dotted curve) of the neutral Higgs bosons as functions of 0i. The values of other parameters 
are the same as in Fig. 2. 

Fig. 3b : The plot of mass matrix elements M15 (solid curve), M25 (dashed curve), and M45 
(dotted curve) of the neutral Higgs bosons as functions of 0i. The values of other parameters 
are the same as in Fig. 2. 

Fig. 4a : The ratios of mass matrix elements |Mi3|/|Mi2| (solid curve), | Afas | / j Af 12 1 (dashed 
curve), and |A^34|/|Mi2| (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 

Fig. 4b : The ratios of mass matrix elements |Afi3|/|Afi4| (solid curve), IM23I/IM14I (dashed 
curve), and IM34I/IM14I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(f>\. The values of other parameters are the same as in Fig. 2. 

Fig. 4c : The ratios of mass matrix elements IM13I/IM24I (solid curve), |M23|/|Ai^24| (dashed 
curve), and |A^34|/|M24| (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 

Fig. 4d : The ratios of mass matrix elements IM13I/IM35I (solid curve), IM23I/IM35I (dashed 
curve), and |M34|/|Af35| (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 

Fig. 4e : The ratios of mass matrix elements |Mi5|/|Mi2| (solid curve), IM25I/IM12I (dashed 
curve), and IM45I/IM12I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 

Fig. 4f : The ratios of mass matrix elements |Mi5|/|Mi4| (solid curve), IM25I/IM14I (dashed 
curve), and |M45|/|Afi4| (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 
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Fig. 4g : The ratios of mass matrix elements IM15I/IM24I (solid curve), IM25I/IM24I (dashed 
curve), and IM45I/IM24I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(pi. The values of other parameters are the same as in Fig. 2. 

Fig. 4h : The ratios of mass matrix elements IM15I/IM35I (solid curve), IM25I/IM35I (dashed 
curve), and IM45I/IM35I (dotted curve) of the neutral Higgs bosons are plotted as a function of 
01. The values of other parameters are the same as in Fig. 2. 

Fig. 5a : The squares of the elements 0^3 (solid curve), O23 (dashed curve), and O34 (dotted 
curve) of the orthogonal matrix defined in the text are plotted as functions of (pi . The values of 
other parameters are the same as in Fig. 2. 

Fig. 5b : The squares of the elements Ofg (solid curve), Ofs (dashed curve), and O45 (dotted 
curve) of the orthogonal matrix defined in the text are plotted as functions of (pi. The values of 
other parameters are the same as in Fig. 2. 

Fig. 6 : The plot of p as a function of (pi. The values of other parameters are the same as in 
Fig. 2. 

Fig. 7 : We plot /Jmax as a function of (pi for A = 300 GeV. The maximum value of p is 
obtained by randomly searching 20,000 points for given (pi in the region of the parameter space 
whose boundaries are set as < (pt = (pb = <Pt = <Pc < 27r, 2 < tan/3 < 40, < A, A; < 0.7, 
< Ax,Ak,x,M2 < 1000 GeV, < At = = < 1000 GeV, < mg = rriL < 1000 GeV, 
and < niT = tub = the < lOOOGeV. 

Fig. 8 : The 20,000 points of Rhi, randomly evaluated in the same parameter region as Fig. 7, 
are plotted against C3, under the condition of (m/i. — nih,.) > 10 GeV (n = 1 to 4). 
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Fig. 1: The plot of the neutral Higgs boson masses as functions of A for (pi = IOtt/IQ. We 
set tan/3 = 5, A = 0.5, k = 0.4, = 590 GcV, Ak = 5 GeV, x = 178 GeV, tuq = rriL 
= 967 GeV, mr = ms = = 794 GeV, At = = = 907 GeV, M2 = 357 GeV, 
4>t = (t>h = (t>T = 4>c = 107r/19. 



32 



mhi(GeV) 



500 



400 



300 



200 



100 




In 

Fig. 2: The plot of the neutral Higgs boson masses as functions of for A = 300 GeV. The 
remaining parameters have the same values as in Fig. 1. 
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Fig. 3a: The plot of mass matrix elements M13 (solid curve), M23 (dashed curve), and M34 
(dotted curve) of the neutral Higgs bosons as functions of The values of other parameters 
are the same as in Fig. 2. 
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Fig. 3b: The plot of mass matrix elements M15 (solid curve), M25 (dashed curve), and M45 
(dotted curve) of the neutral HIggs bosons as functions of <pi. The values of other parameters 
are the same as in Fig. 2. 
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Fig. 4a: The ratios of mass matrix elements |Mi3|/|Mi2| (solid curve), IM23I/IM12I (dashed 
curve), and IM34I/IM12I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(1)1. The values of other parameters are the same as in Fig. 2. 
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Fig. 4b: The ratios of mass matrix elements |Mi3|/|Mi4| (solid curve), IM23I/IM14I (dashed 
curve), and IM34I/IM14I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(f)!. The values of other parameters are the same as in Fig. 2. 
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Fig. 4c: The ratios of mass matrix elements IM13I/IM24I (solid curve), IM23I/IM24I (dashed 
curve), and IM34I/IM24I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(f)!. The values of other parameters are the same as in Fig. 2. 
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Fig. 4d: The ratios of mass matrix elements IM13I/IM35I (solid curve), IM23I/IM35I (dashed 
curve), and IM34I/IM35I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(f)!. The values of other parameters are the same as in Fig. 2. 
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Fig. 4e: The ratios of mass matrix elements |Mi5|/|Mi2| (solid curve), IM25I/IM12I (dashed 
curve), and IM45I/IM12I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(f)!. The values of other parameters are the same as in Fig. 2. 
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Fig. 4f: The ratios of mass matrix elements |Mi5|/|Mi4| (solid curve), |-/Vf25|/|Mi4| (dashed 
curve), and |M45|/|Mi4| (dotted curve) of the neutral Higgs bosons are plotted as functions of 
01. The values of other parameters are the same as in Fig. 2. 
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Fig. 4g: The ratios of mass matrix elements IM15I/IM24I (solid curve), IM25I/IM24I (dashed 
curve), and IM45I/IM24I (dotted curve) of the neutral Higgs bosons are plotted as functions of 
(1)1. The values of other parameters are the same as in Fig. 2. 
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Fig. 4h: The ratios of mass matrix elements IM15I/IM35I (solid curve), IM25I/I-M35I (dashed 
curve), and IM45I/IM35I (dotted curve) of the neutral Higgs bosons are plotted as a function of 
(f)!. The values of other parameters are the same as in Fig. 2. 
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Fig. 5a: The squares of the elements Of^ (sohd curve), (dashed curve), and (dotted 
curve) of the orthogonal matrix defined in the text are plotted as functions of ^1. The values of 
other parameters are the same as in Fig. 2. 



44 



^(Oiil 




Fig. 5b: The squares of the elements Of^ (sohd curve), (dashed curve), and OI5 (dotted 
curve) of the orthogonal matrix defined in the text are plotted as functions of (f)i. The values of 
other parameters are the same as in Fig. 2. 
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Fig. 7: Wc plot Pmax as a function of for A = 300 GcV. The maximum value of p is 
obtained by randomly searching 20,000 points for given (pi in the region of the parameter space 
whose boundaries are set as < (f)t = (f)b = (f)T = (f>c < 27r, 2 < tan/3 < 40, < A, A; < 0.7, 
< Ax,Ak,x,M2 < 1000 GeV, < At = = < 1000 GeV, < mg = rriL < 1000 GeV, 
and < rriT = tub = tue < lOOOGeV. 
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Fig. 8: The 20,000 points of Rhi, randomly evaluated in the same parameter region as Fig. 7, 
are plotted against C3, under the condition of imh^^^_^^^ — fn^^^^) > 10 GeV (n = 1 to 4). 
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